function res = incoming_wave(kh,theta,x,y)

kx = kh * cos(theta);
ky = kh * sin(theta);

i = sqrt(-1);

res = exp(i*(kx.*x + ky.*y));

return